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The Peierls distortions in a two-dimensional electron-lattice system described by a 
Su-Schrieffer-Heeger type model extended to two-dimensions are numerically studied 
for a square lattice. The electronic band is just half-filled and the nesting vector is 
(n/a, it /a) with a the lattice constant. In contrast to the previous understanding on 
the Peierls transition in two dimensions, the distortions which are determined so as to 
minimize the total energy of the system involve not only the Fourier component with the 
nesting wave vector but also many other components with wave vectors parallel to the 
nesting vector. It is found that such unusual distortions contribute to the formation 
of gap in the electronic energy spectrum by indirectly (in the sense of second order 
perturbation) connecting two states having wave vectors differing by the nesting vector 
from each other. Analyses for different system sizes and for different electron-lattice 
coupling constants indicate that the existence of such distortions is not a numerical 
artifact. It is shown that the gap of the electronic energy spectrum is finite everywhere 
over the Fermi surface. 

KEYWORDS: Peierls transition, two-dimensional electron-lattice systems, nesting, Peierls gap, Peierls 
distortion, second order perturbation 

§1. Introduction 

The Peierls transition is caused by the freezing of a lattice distortion mode which can 
connect degenerate electronic states at the Fermi level.0* The presence of such a distortion 
induces an energy gap at the Fermi level of the electronic spectrum. This gap which is 
called Peierls gap lowers the electronic energy. In some cases, this reduction of energy 
overcomes the increase of the lattice energy due to the frozen mode. In one-dimensional 
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systems, particularly, the lowering of the electronic energy is proportional to the square 
times logarithm of the frozen mode amplitude and therefore can overcome the lattice energy 
increase which is proportional to the square of the amplitude. Because of the competition 
between the decrease of the electronic energy and the increase of the lattice energy, there 
exists a value of the frozen mode amplitude minimizing the total energy of the electron-lattice 
system. 

In dimensions higher than one, the situation is not so simple. This is because the number of 
states at the Fermi level is only two in the case of a one- dimensional system, whereas those in 
two- or three-dimensional systems form equi-energy line or surface, respectively. In general, a 
single mode of the lattice distortion can connect only two points in the Fermi surface (or line). 
In this situation, the gain in the electronic energy is too small to overcome the increase of 
the lattice energy. However, in some special situation where a single lattice distortion mode 
can connect many states at the Fermi level, the electronic energy is lowered substantially 
and the Peierls transition becomes possible. This situation is known as "nesting". The 
simplest case can be seen in the two-dimensional square lattice tight-binding model with a 
half-filled electronic band. In this case, the Fermi line is a square within the first Brillouin 
zone combining four points, (ir/a,0), (0,ir/a), (— 7r/a,0) and (0,— 7r/a), with a the lattice 
constant. The nesting vector of this system is Q = (7r/a,7r/a); Q = (rr/a, — n/a) is also a 
nesting vector, but it is equivalent to Q. 

The Peierls distortion was extensively studied in the late eighties to the early nineties in 
connection to the high T c superconductors.i'H , 0'[l A typical argument was given, e.g., by 
Tang and Hirschl who used a tight-binding model where the transfer integral is modified 
by the lattice displacement; its one-dimensional version is known as Su-Schrieffer-Heeger's 
(SSH) model.!) According to the analysis given in ref. ||], the lattice distortion minimizing 
the total energy has the wave vector Q and a polarization parallel to one of the main crystal 
£1X6 S clS shown in Fig. 1. 

As will be discussed in this paper, the straight forward numerical study of the ground 
state of the tight-binding square lattice electron-lattice system described by the 2D version 
of the SSH model indicates that the distortion pattern shown in Fig. 1 does not give the 
lowest energy state. Fourier analysis of the distortion pattern yielding the lowest energy 
state shows the freezing of many other modes with wave vectors parallel to Q in addition 
to that of the main mode with Q. We discuss in this paper why this is the case, along with 
the results of numerical calculations. 

In the next section, the model and the method of determining the lowest energy state are 
described. In § 3, the reason why many modes with wave numbers different from Q can 
contribute to lower the total energy will be discussed. Based on the discussion in § 3, we 
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Fig. 1. The asymmetric dimcrization pattern. The bonds with thick horizontal lines represent shorter 
bonds, and the other horizontal bonds are longer ones. There is no dimerization in the perpendicular 
direction in this situation. 

will propose in § 4 a new method to find the lowest energy state, which can reduce a single 
2D problem essentially to many ID problems. The last section is devoted to summary and 
discussion. 

§2. Model and Basic Formulation 

The model Hamiltonian treated in this paper is the following SSH-type Hamiltonian ex- 
tended to two dimensions, 

H = - J2i(to - ax i:j )(cl +1JtS c iJ:S + h.c.) 

i,j,s 

+(t - aVij)(c\j+i,s c ij,s + h - c -)] 

+f E«+*a (2.i) 

where represents a lattice point in a square lattice with a lattice constant a, and 

j) - u x (i,j), y itj = u y (i,j + 1) - u y (i,j) with u(i,j) the displacement vector 
of an ion-unit at the site to the transfer integral of the equidistant (undeformed) 

lattice, a the electron-lattice coupling constant, c\j s and q JiS the creation and annihilation 
operators of an electron at the site and with spin s. The last term on the right hand 
side of eq. (2.1) describes the lattice harmonic potential energy with K the force constant. 
The lattice kinetic energy is omitted since we do not discuss the dynamical problem in this 
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paper. 

If we assume the Peierls distortion having the nesting vector Q = (Q x , Q y ) = (n/a, it /a) 
as done by Tang and Hirsch,& the bond variables Xij and yij can be expressed in the 
following form, 

Xij = Xoe i (Q* i +Qv3) a = (-iy+i Xo , (2.2) 
Vij = jfoe i(0 " <+Q » J ' )o = {-l) i+j y , (2.3) 

where xq and yo are the amplitudes of the distortion to be determined to minimize the total 
energy of the system. In this situation the electronic part of the above Hamiltonian can be 
easily diagonalized by introducing the Fourier expansion of the electronic field operators as 
follows, 

^• = ^E^e i(fc - + ^>, (2.4) 
iV k 

where N 2 is the total number of lattice points and the components of the wave vector 
k = {k x , k y ) are given by integer times 27r/Na, respectively, on the assumption of the periodic 
boundary conditions. The diagonalization process is straightforward, and the ground state 
energy of the system is expressed in the form, 

E? ! = -2j2 , E k + N^( X l + yl), (2.5) 
k A 

where the sum over the wave vector k is restricted to the region satisfying the following two 
conditions, 

-- <k x + k v < - (2.6) 
a a 

~- <k X - ky< ~ (2.7) 

a a 

and Efc is given by 



E k = V4 + A l (2-8) 
with Efc = —2to(cos k x a + cosk y a) and = 2a(xo sin k x a + yosmk y a). The factor 2 in 
front of the fc-sum in eq. ( |2.5| ) is due to the spin degeneracy. The boundary determined by 
the conditions eqs. (2.6) and (2.7) is nothing but the Fermi surface (or line since the present 
system is two dimensional) and on this boundary vanishes. Across this boundary the 
electronic energy band is divided into two pieces, the dispersions of the lower and upper 
bands being expressed by — and E^, respectively. This means that the energy gap at the 
Fermi surface is given by 2|A^|. 

According to Tang and Hirsch,! 1 the dimerization pattern shown in Fig. 1 which is realized 
by setting x =f and yo = or vice versa and therefore highly asymmetric in x and y 
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directions can have lower energy than the symmetric dimerization realized when x = y /= 
(or equivalently x Q = —yo)- It is n °t difficult to understand the reason if we see the expression 
of Afc for each case; in the asymmetric case, A^ ym = 2ax sin k x a, and in the symmetric case, 
A^ m = 4ax sin[|(fca; + k y )a\ cos[\{k x — k y )a\. The gap on the Fermi surface does not vanish 
except for the points (0, n/a) and {it /a, 0) in the case of the asymmetric dimerization. On the 
other hand the gap vanishes on the line k x — k y = it /a in the case of symmetric dimerization. 
We have confirmed numerically that if we fix the ratio r = yo/x (0 < r < 1) and minimize 
the total energy with respect to Xo, then the minimum value of the energy is a monotonically 
increasing function of r and becomes the smallest at r = 0. 

It is clear that the dimerization pattern shown in Fig. 1 can yield the lowest energy state 
as far as we consider only the distortion with the basic nesting vector Q. However, we should 
note that even in this asymmetric dimerization the gap vanishes at some special points on 
the Fermi surface. In general the Peierls gap is proportional to the matrix element of the 
electron-lattice coupling term in H between two electronic states with wave vectors k and 
k ± Q. In the case of the SSH-type Hamiltonian as used in this paper, this matrix element 
is given by a linear combination of sin k x a and sin k y a. It vanishes at (ir/a,0) and (0, 7r/a) 
irrespectively of the lattice dimerization pattern. This means that the degeneracy between 
(ir/a, 0) and (0, — n/a) [or between (0, ir/a) and (— ir/a, 0)] cannot be removed within the 
first order perturbation due to the lattice distortion with Q. We shall come back to this 
problem in § 4. 

§3. Numerical Study 

In order to check whether the asymmetric dimerization pattern shown in Fig. 1 can yield 
the lowest energy state or whether there exist any different dimerization patterns giving 
still lower energy, we have studied numerically the lowest energy state of the Hamiltonian 
eq. (2.1). Once we know the local values of Xj^-'s and yij's, the electronic wave functions 
{4> u {i,j)} are calculated along with corresponding eigenenergies {e u } from the following 
Schrodinger equation, 

e v <j> v (i,j) = -{t - axij)4) v {i + 1, j) 
-(t - axi_ij)<i) u {i - l,j) 
-{t - ay itj )<i> v (i,j + 1) 

-{to - otyij-i)(f> v (i,j - 1)]. (3.1) 

Since the number of electrons is fixed, it is straightforward to obtain the electronic ground 
state energy for the given configurations of x^s and y^j's. The bond length variables x^/s 
and yi,/s are also involved in the lattice potential energy, and they are determined so as to 
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minimize the total energy of the system. This condition yields the following self-consistent 
equations for these variables similarly as in one-dimensional casesJ3& 

+j^EE^ + 1 -jM.i). ( 3 - 2 ) 

+^E&(i./ + lMj"). ( 3 - 3 ) 

jl V 

where the summation over the one particle states v is restricted to the occupied ones in the 
electronic ground state; note that the spin degeneracy factor should be included. The second 
terms on the right hand sides are due to the periodic boundary conditions, which require 

N N 

Xij = for arbitrary j and ^ = for arbitrary i. 

In most of the practical calculations we use typically the following values of parameters; 
to = 2.5eV, K = 21.0eV/A 2 and a = 4.0eV/A. These values are near to those for poly- 
acetylene. As is well known what is important in this type of electron-lattice systems is the 
dimensionless coupling constant defined by A = a 2 /Kt^. Aforementioned values of param- 
eters give A ~ 0.30. When we want to study the coupling constant dependence of various 
properties, only the value of a is changed for simplicity. 

The set of self-consistent equations (3.1) to (3.3) can be solved numerically by iteration. 
First we give initial values for Xij's and yi,/s, and then calculate (j) u (i,jys by a matrix 
diagonalization subroutine, the results of which are substituted into the right hand sides of 
eqs. (3.2) and (3.3). The resulting values of x^/s and ytj's are compared with the previous 
values. If the difference is not small, we proceed by replacing the initial values of a^j's and 
Vi,/ 8 by new ones. This procedure is continued until the difference becomes negligibly small. 

We try three different initial configurations. First one is the asymmetrically dimerized 
pattern as shown in Fig. 1. If we start the iteration with the uniform dimerization as 
expressed by eqs. (2.2) and (2.3) by setting y = and giving an appropriate value of the 
order of 10 _2 a for xq, we end up with the same pattern with a value of xq which minimizes 
the total energy of the system within that pattern. The second choice is the symmetrically 
dimerized pattern as given by eqs (2.2) and (2.3) with xq = yo ^ 0. The third one is a 
random distortion; we give random values for the lattice displacement vectors {u(i,j)} with 
a relatively small amplitude of the order of 10~ 2 x a; many samples are studied in this 
treatment. The last two choices yield essentially the same result. The resulting distortion 
pattern is not described by eqs. (2.2) and (2.3). An example of Fourier analysis of this 
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pattern is shown in Fig. 2, where X(q x , q y ) and Y(q x , q y ) are defined by 

1 

1 

We have found that all the Fourier components with q x ^ q y vanish® 



X(q x ,q y ) = -^Y. x id exp[-i(q x i + q y j)a], (3.4) 
y{Qx,q y ) = -M^yijexpl-ifaJ + qyfia]. (3.5) 
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Fig. 2. The Fourier spectra of the distortion pattern obtained by numerical minimization of the total 
energy. As for the direction of the wave number q, the Fourier components vanish when q x =/= q y . The 
abscissa indicates q^a (— q x a — q y a). The system size is 36 x 36. The electron-lattice coupling constant a 
is equal to 4.0eV/A 2 which corresponds to A = 0.30. 



What is characteristic in these Fourier spectra is that the amplitudes [= 
\J\X(q^)\ 2 + |y(qj)| 2 {i = 1,2)] of two modes with q 1 and q 2 satisfying the condition 
Qi + <?2 = Q are equal to each other. The meaning of this property will be discussed 
in the next section. It should be noted that as far as the Q-components are concerned the 
amplitudes are highly isotropic in x and y directions. 

The electronic energy spectrum corresponding to this lattice configuration is different from 
that in the asymmetric dimerization case (Fig. 1). There is no point on the Fermi surface 
where the Peierls gap vanishes. The behavior of the gap along the Fermi surface will be 
discussed later in detail. 
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In order to check whether this complicated state has a lower energy than the asymmetri- 
cally dimerized state as shown in Fig. 1, we study the size dependence of the energy difference 
between the asymmetrically dimerized state (Fig. I) and the lowest energy state obtained 
by solving the self-consistent equations (3.1) to (3.3); the total energy of the former state is 
denoted by Ea and the latter by E G) both of them being negative. The result is summarized 
in Fig. 3, where the energy difference scaled by |£a| is plotted as a function of N~ 2 , the in- 
verse of the total number of the lattice points, for two different values of a, other parameters 
being fixed as mentioned before. 




Fig. 3. The size-dcpcndcnce of the energy difference between the asymmetrically dimerized state and the 
lowest energy state obtained by solving the self-consistent equations. The abscissa is the inverse of the total 
number of the lattice points. The value of a is indicated in the figure along with the corresponding value 
of the dimensionless coupling constant A, other parameters being fixed (to = 2.5eV and K = 21.0eV/A 2 ). 
The continuous lines are fitting to quadratic polynomials. 



From Fig. 3, we may safely conclude that the energy difference remains finite in the ther- 
modynamic limit, and that the lowest energy state obtained from the self-consistent equation 
has certainly a smaller energy than the asymmetrically dimerized state shown in Fig. 1. 

§4. Peierls Gap due to Second Order Process 

In this section we discuss why the single mode dimerization (Fig. 1) is not the lowest 
energy state. The lowest energy state obtained numerically in the previous section involves 
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multi-mode distortions. The Fourier analysis of the distortions indicates that the gap in 
the electronic spectrum is induced not only by the first order perturbation but also by the 
second order process. As discussed in § 1, the Peierls gap in the electronic spectrum is 
formed usually by the first order process in which a state \k) on the Fermi surface is coupled 
to another state \k ± Q) on the Fermi surface by the lattice distortion with a wave vector Q 
(the nesting vector). The consideration developed in § 2 indicates that the gap vanishes at 
special points on the Fermi surface, (ir/a, 0) and (0, vr/a), by this first order process because 
of the peculiar wave number dependence of the electron-lattice coupling term. This fact 
means that in order to get a finite gap at those point we have to consider the second order 
process where at least two lattice distortion modes should be relevant. 

If two distortion modes with wave vectors q 1 and q 2 satisfying q 1 + q 2 = Q are involved in 
the second order process, there should appear matrix elements connecting one of the states 
|fe) or \k + Q) and one of \k + q ± ) or \k + q 2 ), where the former states are on the Fermi 
surface. This indicates in turn that the two states {k + q^ and \k + q 2 ) which are not on the 
Fermi surface are mixed with the states |fc) and \k + Q) on the Fermi surface and contribute 
to the formation of the gap at the Fermi surface. In such a situation it is natural to expect 
the creation of lattice distortions with wave vectors q l — q 2 , q 2 — q 1 , q 1 — q 2 + Q{= 2q 1 ) 
and q 2 — q 1 + Q{= 2qr 2 ) . In this way, many lattice distortion modes can be involved in the 
second order process. 

Based on the numerical results discussed in the previous section, we assume that, in the 
lowest energy state, only the lattice distortions with wave vectors parallel to Q are existing. 
Namely the bond variables Xij and y^j are assumed to be expressed in the form, 



By this assumption the lattice degree of freedom is reduced from N 2 to iV — 1; the uniform 
mode with q = is excluded trivially. Then the original Hamiltonian can be written in the 
wave number representation as follows, 




(4.1) 



0<q<n/a 



y w = ItoC-l)*'" + E [y^^+c.c.]. 



(4.2) 



0<q<ir/a 



H = H £ k%s c K 



k,s 



,s 




k,s 




0<q<n/a k,s 
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+y*cos (fcy-|) a] 4_ g s c fe 

+iV 2 f(^ + l/o 2 ) 

+N 2 K £ (W + 1%1 2 ), 



} 



(4.3) 



0<g<7r/a 



where the vector q stands for (g, g). 

Thus, in order to obtain the electronic energy spectrum in the presence of the lattice 
distortions considered above, we have to solve the eigenvalue problem of the following form, 



where ip(k x , k y ) means the wave function in the wave number representation. It should 
be noted that, in the above equation, only the combinations of k x and k y with a constant 
difference are involved. It is not difficult to understand from this fact that the Hamiltonian 
matrix whose original size is N 2 x iV 2 can be decomposed into N pieces each of which has 
a size N x N. Therefore the most general eigenvalue problem eq. (3.1) is decoupled into 
N one dimensional problems in the present situation. Each eigenvalue equation of the form 
eq. (4.4) yields a set of N eigenvalues and N eigenfunctions. By solving N sets of eigenvalue 
problems for one set of {x q , y q } and (xq, yo), we end up with a set of N 2 eigenvalues and 
N 2 eigenfunctions. Sorting N 2 eigenvalues in increasing order, we determine the state index 
v and attach it to each eigenvalue and eigenfunction as e v and ip v {k + k', k). If we note the 
periodicity of the wave functions in the fc-space, we have only to consider the range of k and 
k' such as —n/a < k, k' < n/a. 

Next we have to discuss how to determine the set of variables, {x q , y q } and (xq, yo). They 
are determined so as to minimize the total energy of the system as done in § 3. The self- 



eip(k x ,k y ) = e k ip(k x ,k y ) 



— 2ia(x sin k x a + y sin k y a) x 

1p{k X ~ Qxi ky — Qy) 




(4.4) 
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consistent equations for these variables are written in the following form, 

' KN 2 - 



x = — JZ' sin k x aipt(k x + 7r/ a, k y + n/a) x 



^A k x,k y ), (4.5) 
y ° = ~ Kiv 2 ^ ^ sin k y a ^ kx + 71 1 a ' ^ + ^/a) x 

" k 

4>u(k x ,k y ), (4.6) 



— > > cos Kj. ax 

i^V 2 V V 2 ' 



ipt(k x - q,k y - q)ip u (k x ,k y ), (4.7) 



2e i<?a/2 



a Y!Y, cos (ky-7; 
" k V Zy 

Vv{kx~ q,k y - q)i/j„(k x ,k y ), (4.8) 



v k 



where the sum over v is restricted to the occupied states. We have omitted complex conjugate 
forms of eqs. (4.7) and (4.8). It should be noted that for each v the corresponding wave 
function is finite only when k x — k y =const., though this constant depends on v. 

In the present calculation scheme, we have only to diagonalize iV different N x N matrices 
instead of a single iV 2 x N 2 . This makes the computational load much lighter, and therefore 
we can treat far larger system than in the case where we deal with two dimensional problems 
directly. 

In Figs. 4 and 5 we show the wave number dependences of the amplitudes and the phases 
of Xg's and y 9 's which have been obtained by solving the self-consistent equations for the 
system size iV = 128 (the total number of lattice points = 128 x 128) and the coupling 
constant a = 4.0 eV/A (A = 0.30). The initial condition used in this iterative calculation is 
Re x q = Im x q = x- m (some constant of the order of 10~ 2 a) for q = 2-rr/Na, x q = for other q 
and y q — for all q, where q is restricted in the region < q < it /a, and xq = yo = x- m . The 
final results does not depend on the values of X{ n . The results for different types of initial 
conditions will be discussed later. 

The amplitudes of x q and y q are completely equal to each other and are found to vanish 
at q values which are even integer times 27r/Na, though, as a matter of course, q = n/a is 
exceptional. As for the phases, we find the relation arg(a;g 1 )+arg(a;g 2 )=arg(?/ (?1 )-|-arg(t/g 2 )=0 
mod(7r) for qi and q^ satisfying q 1 + q 2 — ir/a. This will be reasonable if we remind us the 
second order perturbation mechanism of the Peierls gap formation. 

Once we know the values of {x q }, {y q }, Xq and yo, we can obtain N eigenenergies for each 
set of N electronic wave vectors. Assuming that the energy of the lower (higher) band is an 
increasing (decreasing) function of the distance from the line k x + k y = 0, we can assign the 
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Fig. 4. The q dependences of the amplitudes \x q \ and \y q \. The values of \xq\ and \yo\ are also plotted. The 
system size is N = 128, and the coupling constant is a = 4.0eV/A (A = 0.30). 



energy versus k relation. The dispersion relation obtained in this way is shown in Fig. 6 for 
the case with a = 6.0eV/A (A = O.69).0> Only the lower band which is fully occupied in the 
ground state is shown along the lines depicted in the sub-figure. Because of the electron-hole 
symmetry of the system the dispersions of the unoccupied levels are the same as those shown 
in Fig. 6 except for the sign of energy. 

Thus by assigning the dispersion relation for all the points in the fc-space, we can now see 
the gap structure on the Fermi surface. As will be clear from Fig. 6, the gap is constant 
along the line (— n/a,0) — (0, 7r/a). On the other hand the gap structure along the line 
(0, n/a) — (n/a, 0) is not so simple as will be found from Fig. 7, where the dispersion curves 
of upper and lower bands along the line (0,n/a) — (it/a,0) are shown for the cases with a 
couple of different values of the electron-lattice coupling constant. 

In contrast to the asymmetrically dimerized case (Fig 1), the gap does not vanish at any 
point, taking almost the same value as that at (tt /2a, it /2a). It will be noteworthy that the 
gap is not minimum at k — (0, 7r/a) or (7r/a, 0) and that the multi-mode effect is almost 
negligible around the point k = (n/2a, 7r/2a). The position of the gap minimum depends on 
the coupling constant. In the same figure, the dispersions in the asymmetrically dimerezed 
case for each coupling constant are shown for comparison. 

In order to see the size dependence of the gap A, we plot in Fig. 8 the gap at three different 
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Fig. 5. The q dependences of the phases aig(x q ) and arg(y 9 ). The values of arg(a;o) and arg(y ), which are 
zero in this example, are also given for comparison. The system size is N = 128, and the coupling constant 
is a = 4.0eV/A (A = 0.30). 



points on the line (0, it /a) — (ir/a, 0), as indicated in the inset, for a coupling constant A = 0.69 
(a = 6.0eV/A). From this figure we find that N = 128 belongs to the large size limit as far 
as the behavior of the gap concerns and that in the large size limit the gaps Ai and A2 at 
(n/2a, 7i/2a) and (0, n/a) are indistinguishable though the value of A 3 remains smaller than 
A x and A 2 . For smaller values of the coupling constant, a similar tendency is found, but 
the large size limit is seen only at larger system sizes, although N = 128 is sufficiently large 
even for a = 4.0eV/A. 

As indicated by the results shown above, the lowest energy state obtained here is not 
completely symmetric with respect to x and y directions. This fact means that there exist 
at least two degenerate states where the roles of x ant y are interchanged. 

In order to check other degeneracy of the lowest energy states, we have studied what kind of 
distortion patterns could be obtained when we change the initial distortions in the iterative 
calculation. As a result, we got many different distortion patterns with the same energy as 
that of the distortion pattern shown in Figs. 4 and 5. Some are different only in the phases of 
the Fourier components of the distortion. Some show completely different behaviors. Among 
them there is a pattern where the nonvanishing Fourier components are x 7T / 2a , 2/vr/2a, ^0 an d 
yo and other components are zero. In these degenerate states, not only the total energy but 
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Fig. 6. The dispersion relations of occupied levels for a = 6.0eV/A (A = 0.69) and the system size 128 x 128. 
The wave numbers are changed along the lines indicated in the sub-figure. 

also the electronic part and the lattice potential part are also the same, respectively. At the 
moment we cannot say how many states are degenerate. Detailed analysis of the degeneracy 
problem is left for a future work. 

§5. Summary and Discussion 

The lowest energy state of a two dimensional electron-lattice system with a half-filled 
electronic band and with a square lattice structure is studied within the SSH-type model 
extended to two dimensions. On the contrary to the previous common understanding that 
only the lattice distortion with the nesting vector Q = (ir/a, tt/o) is frozen in the lowest 
energy state,!* many modes are found to be frozen in the real lowest energy state. The 
state discussed by Tang and HirscfM might be a local minimum but it is not the absolute 
minimum of energy. This is because the Peierls gap vanishes at (iv/a, 0) and (0, 7r/a) due 
to the wave number dependence of the electron-lattice coupling term. We have pointed out 
that the second order perturbation mechanism of the Peierls gap formation is important. 
Numerical minimization of the total energy leads to the conclusion that many modes having 
the wave number parallel to Q contribute to the formation of the gap. Assuming this is the 
case, we can reduce the two-dimensional problem into one-dimensional problems by using the 
wave number representation of the Hamiltonian as discussed in § 4. In this formulation it is 
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Fig. 7. The dispersions of the highest occupied and lowest unoccupied levels along the line (0, it /a) — 
(w/a,0). Three different values of electron-lattice constant are considered, a[eV/A] ~ 4.0 (A = 0.30); 
continuous thick line, 4.4 (0.37); dotted line and 5.2 (0.52); dash-dotted line. The similar dispersions for 
the asymmetric dimerized state (Fig. 1) are also shown by the same type thinner lines. 



possible to treat far larger system sizes than in the case where we search a minimum energy 
state of the two-dimensional system directly. The electronic energy dispersion and the gap 
structure on the Fermi surface have been analyzed. Although the wave number dependences 
of the amplitudes and phases of condensed modes look to show a certain symmetry in the x 
and y directions, the electronic structures are not necessarily symmetric. 

What we have shown in this paper is not the unique lowest energy state. A preliminary 
study indicates there might be infinite number of degenerate ground states. The degree 
of degeneracy is not known at the moment. It is not clear also what kind of symmetry is 
relevant to this degeneracy. Detailed study of the degeneracy problem is left for the future 
work. Nevertheless it will be worthwhile to mention that this type of degeneracy of the 
ground state yields the possibility of the formation of domain walls like the solitons in the 
one-dimensional systems where the number of degenerate ground state is only two. 

100) 

This kind of domain walls connecting two different ground states may supply interesting 
properties of the two-dimensional electron-lattice systems just as the charged and neutral 
solitons in polyacetylene did. 



15 




Fig. 8. The system size dependence of the gap A in the electronic spectrum at different points on the line 
(0,n/a) — (n/a,0). The coupling constant is fixed at A = 0.69. 
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